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Abstract 

We derive the effects of a non-zero cosmological constant A on gravitational wave 
propagation in the linearized approximation of general relativity. In this approximation 
we consider the situation where the metric can be written as g^jy = 77^1/ + h^^ + /ij^, 
h^i)^ << 1, where h^^ is the background perturbation and /ij^ is a modification inter- 
pretable as a gravitational wave. For A 7^ this linearization of Einstein equations is self- 
consistent only in certain coordinate systems. The cosmological Friedmann-Robertson- 
Walker coordinates do not belong to this class and the derived linearized solutions have 
to be reinterpreted in a coordinate system that is homogeneous and isotropic to make 
contact with observations. Plane waves in the linear theory acquire modifications of order 
\/A, both in the amplitude and the phase, when considered in FRW coordinates. In the 
linearization process for we have also included terms of order O(Ah^iy). For the 
background perturbation h^^^ the difference is very small but when the term h^^^A is re- 
tained the equations of motion can be interpreted as describing massive spin-2 particles. 
However, the extra degrees of freedom can be approximately gauged away, coupling to 
matter sources with a strength proportional to the cosmological constant itself. Finally 
we discuss the viability of detecting the modifications caused by the cosmological con- 
stant on the amplitude and phase of gravitational waves. In some cases the distortion 
with respect to gravitational waves propagating in Minkowski space-time is considerable. 
The effect of A could have a detectable impact on pulsar timing arrays. 
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1 Introduction 



The smallness of the cosmological constant obtained from fits to the current ACDM cosmolog- 
ical models^l] (A ~ 10~^^ m~^) may lead us to believe that it is totally unobservable except 
at the largest distances. However, the issue of the relevance of the cosmological constant in 
local measurements (meaning measurements that involve sub-cosmological scales, such as for 
instance galaxy clusters) has received growing attention |2l E]. One interesting possibility is 
assessing the influence of A on the bending of light from distant objects. At present there are 
rather diverging results on the subject giving rather different results concerning the relevance 
of A ranging from zero[4j or very small [5] to appreciable ones [6]. The effect of A on the 
photon propagation, including frequency shift, Shapiro time delay and deflection of light, is 
currently under consideration |7j. 

The importance of these studies cannot be overemphasized. The presence of a non- 
zero cosmological constant contributing around 70% to the energy and matter budget of 
the universe, seemingly making the Universe globally a de Sitter space-time, is one of the 
intriguing puzzles of Physics in our time. Observations capable of confirming or refuting the 
relevance of A at redshift z < 1 are clearly of utmost importance. 

The studies of what has been termed 'local gravity with a cosmological constant' rely on an 
approximate solution, valid at first order in A, obtained after linearizing Einstein equations. 
These solutions have recently been studied in detail by one of the authors [3j using different 
gauge choices. It has been found that in the Lorenz gauge one can in addition require time 
independence of the metric solutions. After an additional coordinate transformation these 
solutions correspond to the linearized version of the Schwarzschild-de Sitter exact solution 
of Einstein equations. The modification to the Newtonian limit in such coordinates was also 
discussed in detail [3J. There are some subtleties related to the physical interpretation of the 
different coordinate systems that we shall review below. 

Here we propose to study a different problem. Namely, how A influences the properties of 
gravitational waves (GW). As of today, gravitational waves are an unambiguous prediction of 
General Relativity that has not been tested directly. They are 'observed' indirectly as they 
are the missing ingredient needed to restore the energy balance of some astrophysical binary 
systems [S]. There are three types of experiments potentially capable of yielding a non-zero 
signal in the coming years. Let us summarize their physical and astrophysical reach here: 

Ground based GW detectors such as LIGO [9] can reach sensitivities down to ~ 10"^'^ 
with optimal sensitivity in the region between 10 Hz and 10'^ Hz. The space mission LISAp^ 
will reach a similar sensitivity in the range 10~^ Hz to 10"'^ Hz but will actually be able to set 
relevant bounds on a more extended range of frequencies. Finally the International Pulsar 
Timing Array project[llj or the Square Kilometer Array project[12] are sensitive to lower 
frequencies u < 10~^ Hz but reach only a sensitivity of ~ 10^^'^ going up to ~ IQ-^^ for 
u ~ 10"^^ Hz. These sensitivity ranges are targeted to specific astrophysical phenomena and 
are expected to provide detectable signals and confirm the existence of GW in the coming 
decades. 

Given the present difficulties in asserting the very existence of GW it may seem academic 
to try to find modifications due to the presence of a cosmological constant that is small. 
However, it should be borne in mind that in the inflationary epoch the value of A was much 
larger than at present so these effects might be of relevance for primordial GW. As we will 
discuss in this work the effect of A could be of some relevance for GW travelling very long 
distances and for pulsar timing array projects. On the other hand, some of the results 
presented here we believe are of interest to understand the issue of the gauge choice in the 
presence of A for the linear theory. Finally, it seems interesting in its own right to attempt to 
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understand wave propagation in de Sitter space-time if A is indeed a fundamental parameter 
of nature. 

This paper is organized as follows. In section 2 we discuss the linearization of Einstein 
equations, including a discussion on different gauges and how they affect the wave equation 
for the gravitational field /i^jy. In section 3 we discuss different coordinate realizations of de 
Sitter space-time and their relation. In section 4 we construct background solutions retaining 
terms of order Kh^y. This discussion is extended in section 5 to include GW solutions that 
'feel' the presence of A. In section 6 we analyze the detectability of the effects previously 
calculated. In section 7 we summarize the conclusions of this study. 

Some of the subjects discussed here appear to have received little attention in the past 
although there is an extensive literature on gravitational waves [13] . The effect of A on GW 
has been considered in [141 15]. Physical consequences appear to have been extracted in 
the context of primordial gravitational waves [16j and only indirectly in what concerns the 
evolution of the modes and the power spectrum. 



2 Linearization in the presence of A 

Einstein equations, derived from the Einstein-Hilbert action, read 

Rtiv - ^9fiuR + ^gfiu = -KTf,u (1) 

where R^,y is the Ricci tensor for g^,^ , A > is the cosmological constant and kT^i/ is the source 
term. T^y is the usual stress-energy tensor of matter in the gravitational field generated by 
Qfj^u and K is the dimensionful constant coupling matter and gravity. However, throughout this 
work we will consider T^i, = unless otherwise specified. The inclusion of the cosmological 
constant term leads to curvature even in the absence of any source 

R = 4A. (2) 

We consider the linearized theory where the metric is written as 

5/^1/ — ^iMu + f^^i/j (3) 

rjfj^i, being the Minkowski metric and h^i, << 1. The Ricci tensor to first order in the small 
perturbation h^^, reads 

R^, = ^ (ah^, + h^, - hl^^^ - hl^^^ , (4) 

indices being lowered and raised with r/^j^ and h = rj'^'^hf^u. The theory is invariant under 
coordinate transformations ^ x^^ + (,^{x). For infinitesimal transformations the pertur- 
bation metric h^i, transforms as h^i, — )• h'^^, = h^y + dfj_£^y + dy^^. A gauge choice is possible, 
amounting to selecting a particular class of coordinates, and in fact such a choice is necessary 
if the perturbation h^j^y is to be quantized. In order to discuss GW two different gauge choices 
are particularly appropriate. 

2.1 Lorenz gauge 

In order to describe perturbations around flat space-time it is customary to employ the Lorenz 
gauge. 

d^K = \duh, (5) 
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or 

d^ht^ = 0, (6) 

where ^ 

h^u = h^y - -rjfiuh (7) 

is the trace reversed version of h^^. 

In this gauge, expression is simphfied 



Rfiu — i^^h^y, (8) 

and we obtain the equation of motion 

□ {hfiu - \'n^luh^ + 2A V = -'^^^l^lu (9) 

which has always to be considered together with the Lorenz gauge condition 

Whether the term of order 0{hK) has to be considered or not depends on the relative 
magnitude of h and A. There will be situations where the inclusion of this term is justified 
and may lead to observable consequences. We shall postpone the rest of the discussion on this 
issue to sections 4 and 5. Note, nonetheless, that if the Ah^u term on the l.h.s. is omitted 
(and only in this case) there is a residual gauge freedom within the Lorenz gauge. If we 
perform a linear coordinate transformation 

_^ ^'1^ = + (10) 

equation ([5]) is fulfilled as long as ^'^ is an harmonic function, i.e. □.^'^ = 0. These residual 
coordinate transformations are sometimes termed 'coordinate waves' for rather obvious rea- 
sons. Note also that whether this is a symmetry of the equations of motion or not, depends 
on the terms retained in the linearization; the term A/i^^ breaks this residual coordinate 
invar iance. 



2.2 A gauge 

It will be useful to consider an alternative gauge choice jl7j. which we will term A-gauge. This 
is given by the gauge condition 

d^K = -Ar/.M^^- (11) 
In this gauge the linearized equations of motion look slightly different 

□ - lv,.uh^ - 2AV = 0. (12) 

In particular we note that the term independent of h^i, on the r.h.s. of Q is absent. There 
is a set of coordinate transformations that can be performed without leaving the gauge orbit 
(jlip ; these are transformations x'^^ = x'^ + with 

□^/^ = -A^'^. (13) 

However, in the A-gauge this residual coordinate transformations are not a symmetry of the 
equations of motion regardless of the terms retained in the linearization and therefore cannot 
be used to remove degrees of freedom. Generally speaking, linearization leaves global Lorentz 
transformations as the only symmetry of the equations of motion. The Lorenz gauge is in a 
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way special as some additional freedom to perform local coordinate transformations remains 
if the term Kh^y is neglected. The situation in the A- gauge, where there is no residual 
symmetry, is, on the contrary, the generic one. 

The connection between the two gauge choices in the linear theory is easily made when 
the terms h.h^,v are omitted. It is implemented via the following change of coordinates 

This change of coordinates transforms a solution of D/i^i/ = in the A-gauge (coordinates 
x) to a solution of D/i^jy = — 2Ar/^,y in Lorenz gauge (coordinates x'). Note the simplicity of 
the equation for linear perturbations in the A-gauge if the term of order Kh^^ is omitted. All 
reference to the cosmological constant is eliminated. 

The previous discussion in the A-gauge reminds us that in general, in the linearized 
approximation, the perturbation metric h^^ is expected to have up to six full degrees of 
freedom. Only in certain cases a residual gauge freedom can be used to further reduce the 
number of degrees of freedom. 

Let us elaborate a bit more on this issue as it is conceptually important. In the Lorenz 
gauge, with the term hh^y omitted, the residual symmetry (fTO]l allows us to move freely 
between different coordinate systems, say x' and y', which are not trivially related by Lorentz 
transformations and yet preserve the form of the equations of motion. On the contrary, if 
we undo transformation (jl4p we get two coordinate systems x and y in which the A-gauge 
condition is fulfilled but at best only one of these gauge transformed coordinate systems obeys 
the linearized equations of motion in the A-gauge; the other one is off-shell. That is to say, 
the number of independent degrees of freedom seems to be larger in the A-gauge. However, 
since this is purely due to a gauge choice, the additional apparent degrees of freedom cannot 
correspond to physical ones. 

If the term of order Kh^y is retained, i.e. in the Lorenz gauge the term 2Kh^y on the 
l.h.s. of dH]) or the analogous —2kh^y in the A-gauge are kept, there is no residual symmetry 
whatsoever. Let us take for example ([9]) in the Lorenz gauge; as we will see in detail in section 
5 this generates a genuine mass term and therefore more physical degrees of freedom appear 
associated to h^y. This is not a gauge artifact. 



3 De Sitter space-time 

De Sitter space-time can be described by many coordinate systems. A convenient choice of 
coordinates is Schwarzschild-de Sitter (SdS). These provide a time-independent metric in a 
gauge that is none of the two previously discussed 

and clearly shows the presence of the de Sitter horizon. We note that this metric admits an 
expansion in integer powers of A. Note also that in this metric the spatial part does not quite 
correspond to spherical coordinates. 

At the opposite extreme, one can select a metric that depends only on time and is position 
independent. It is the Friedmann- Robertson- Walker (FRW) metric 

ds^ = dT^ - e-}ip{2^T)dX^. (16) 



A 2 
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This metric incorporates the physical principles of cosmological homogeneity and isotropy 
as it does not depend on the position. The coordinates X'^ have a clear physical meaning, 
they are comoving coordinates anchored in space that expand with the universe. These are 
the natural coordinates where our world appears homogeneous and isotropic. It is easy to 
see that the FRW metric does not fulfill any linearized Einstein equation, even for very early 
times t « when is very close to the Minkowski metric. In fact, no metric that depends 

only on time can be a solution of the linearized Einstein equations; incompatibilities appear 
immediately for any gauge choice. 

One should therefore accept that the linearized Einstein equations in the presence of A 
cannot be imposed in the physically relevant comoving coordinate system. This of course 
has implications on GW as the very concept of 'wave' does require a wave equation, which is 
just impossible in FRW coordinates. On the other hand, the wave equation Oh^i, = found 
in the A-gauge is expressed in a set of coordinates whose meaning is yet to be interpreted. 
Therefore the simplicity of this equation is deceiving. 

We will argue in the next section that the coordinates implied by the choice of the A- 
gauge or of Lorenz gauge are closely related to SdS coordinates. Then the way to proceed 
is to find a solution for GW in the Lorenz gauge, a coordinate system where linearization of 
the Einstein equations is consistent, and then transform the solution to FRW coordinates in 
order to extract observable consequences. 

Both the SdS metric and the FRW metric are valid (but rather different) descriptions of 
de Sitter geometry. One can work out the exact transformation between the two coordinate 
systems 




V3 



+ T 



(17) 



where T and R are respectively the cosmological time and comoving coordinates whose phys- 
ical realization is clear. This transformation is valid inside the cosmological horizon, i.e. 
R< Applying ([HD to ^ we obtain 



(18) 



Now it is immediate to see that the FRW metric does not fulfill any linearized Einstein 
equation, even if i << 1/\/A as it is not expandable in integer powers of A. The same 
transformations for the linearized version of the metrics gives 



ds^ 



ds^ = dT^ 



1 

1 r 



dt^ 



A 2 
1 + -f2 



A. 



1 + 2W- T + 2-T' 



{dR^ + R^dn'^ 



(19) 



which will only reasonably approximate the expansion of FRW for values of i? ~ T < < . 



Note that, although the last metric in (jl9p is linearized, it does not fulfill any linearized 
Einstein equations. 

The previous transformation provides the relationship between a framework where the 
Einstein equations can be consistently linearized and the actual coordinate system in which 
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we observe. The solutions easily found in the linearized theory have to be transformed to 
the physically meaningful coordinate system in order to make predictions. It is at this point 
that non-trivial effects related to A will appear. They are discussed in section 5. Of course, 
given the current value of A, these effects will be small. We believe nonetheless, that these 
corrections are conceptually important. Note also that ()17p involves \/A and not A, yielding 
corrections that are potentially much more relevant for observation than those of order 0(A). 

Equation (jl6p is just one of the many possible cosmological FRW metrics. Other possibil- 
ities such as a power law cosmological scale factor do not correspond to a de Sitter space-time 
and therefore there is no obvious change of coordinates that allows to reexpress a GW, i.e. a 
solution to a wave equation, in that physically meaningful coordinate system. 



4 Background solutions 

We shall work consistently in the linearized approximation both for the background modifi- 
cation h^j^ and for gravitational wave perturbations /ij^. Namely, the metric can be written 
as g^u = V^iv + h^u + ^ju!-) where h^l!^ <C 1. To keep the notation simple we shall only use 
the superscript A when confusion with wave perturbations /ij^ is possible. In this section we 
will be concerned with background linearized solutions when the cosmological constant A is 
present. 

The value of the cosmological constant has presumably not been the same throughout the 
history of the universe. In early epochs, perhaps following an inflationary period, its value is 
believed to have been much larger [18] . This fact suggests that it may be necessary in some 
circumstances to retain the term Ah^^. Likewise it will be necessary for consistency to keep 
terms of order A/iJ^ as the magnitudes of and A are unrelated. 

In what follows we proceed without making any assumptions on the value of A; we will 
just assume that the perturbation that induces on the background metric /i^^ is small enough 
for the linearized approximation to be meaningful. 



4.1 Lowest order solutions 

First we turn to the lowest order solutions already discussed in [3J, which correspond to 
neglecting terms of 0{Kh^y). In the Lorenz gauge this amounts to solving the following 
equation 

□^A*!/ = - 2Ar/^^ 



d^,K =0. 



(20) 



Linearization limits the validity of the solution to values of the coordinates such that « 
1/A. 

Before discussing the solutions to (pOj) we take a look at the equations in the A-gauge 



r (21) 
d^hi^ = - Ar/^^x^. 

Note once more that the linearized equations are not invariant under gauge transformations. 
In the Lorenz gauge the cosmological constant is regarded as a gravitational source, it ap- 
pears in the equations of motion, whereas in the A-gauge all dependency in the cosmological 
constant at this order appears through the gauge condition only and in a way it can be in- 
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terpreted as a consequence of the coordinate choic^. The connection between the two gauge 
choices in the hnear theory has already been discussed. 

We can easily solve equations ([2T|) to find the traceless solution 



h 



A 
18 



(22) 



If we require that the solution is proportional to A and involves only the coordinates x^^ this 
is the unique solution. In addition, (I22p is the only one that is Lorentz-covariant (note that 
r]^y is the underlying metric and there is no other four- vector at our disposal). 

It is worth noticing that since there is no residual freedom in this gauge, no transformation 
can turn this solution into a static metric: The K-gauge is explicitly incompatible with the 
solutions being static. 

We now transform the solution back to the Lorenz gauge using ()14p . We find 



A 



(23) 



Without the A/i^j^ term the equation of motion is actually invariant under residual trans- 
formations. The number of physical degrees of freedom therefore is reduced to two. This is 
the only covariant-looking solution in the Lorenz gauge but only one of the infinite number 
of solutions reachable by non-covariant residual transformations. The most general form of 
such transformations is 



( A{t^ + r'^)t \ 

{Bit"^ + B2X^ + Bsiy^ + z^)) X 
{Bit^ + B2y^ + Bsix^ + z^)) y 
V {Bit^ + B2z^ + Bsix^ + y^)) z j 



(24) 



where 2B\ — 6B2 — 4i?3 = 0. In particular we find the values of these constants that allow us 
to reproduce the static solution of [3j . 



18' 



Si 



, Bi. = — ; 

18' 36' 



(25) 



One should ask at this point what are these coordinates. We already know that they cannot 
correspond to cosmological coordinates. In fact the resulting metric is neither homogeneous 
nor isotropic although it preserves the symmetry among the three axes. The answer becomes 
obvious once one discovers that one of the possible residual gauge transformations eliminates 
the time dependence of the metric. A generalization of Birkhoff 's theorem [19] states that 
there is a unique static solution with spherical symmetry which is the Schwarzschild-de Sit- 
ter metric previously discussed, or more precisely the first order of it in the A expansion. 
Since Schwarzschild-de Sitter does not fulfill the Lorenz gauge condition, a time- independent 
coordinate transformation must also be involved. Let us explicitly show this point using a 
succession of coordinate transformations linear in A. 

The first step is to transform (I23p to a static solution. We start from 



A , o 
1 + -(3*2 
9^ 



2r' 



2A 



1 



A 
9< 



-2f + 2r2 + x^^l 



dx' 



2A 

1 X* dt dx' + — X* x^ dx* dx^ 

9 9 



(26) 



^ This of course does not mean that the consequences of A can be removed by a wise coordinate transfor- 
mation but it does mean that it disappears from the equations of motion themselves. 
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where i = 1,2,3 and i ^ j- After the following change of coordinates 



X =x + 



y=y + 



A 
9" 
A 



./2 



T ^ 4 
2 + 4 



A 



9 
A 
18 



/2 



+ 2/'2) 



(27) 



the metric transforms into the static solution to order A found in [3], 

A 





\ A ,2' 
1 r 


dt'^- 




CO 





l-:^(r-+3x-) 



/2 



(28) 



Note that this solution is still in the Lorenz gauge; we only performed a residual gauge 
transformation that is allowed in this gauge. Since our starting solution is only valid to order 
A, in any change of coordinates, either exact or linear, we only keep terms linear in the 
cosmological constant. We can further transform (j28p to obtain a fully spherically symmetric 
solution. Under the following change 



//3 



X =X + —X 
I // I A 112, 

y =y + Y^y 

z =z H z 

12 



(29) 



t' =t", 



we obtain 



ds'^ = 


1 r 


dt"^ - 


1 r 




CO 




6 



(30) 



which does not obey (j20p anymore. We can now perform another coordinate transformation 
to obtain the SdS metric to order A 



. A „3 
12 



t" =t 



ds' 



"1 ^-2' 


dp - 


A 


1 r 


1 + -f2 


3 




3 



df2 + f^dJl' 



(31) 
(32) 



This is the linearized Schwarzschild-deSitter metric. Essentially the background solution (j23p 
is the SdS metric in a set of coordinates related to SdS by time independent transformations. 



4.2 Next-order solutions 

Let us now relax the approximation of the previous section and retain terms proportional to 
Ah^iy. In particular we will be interested later in terms of order A/i|^ that will influence the 
propagation of gravitational waves. 

In the Lorenz gauge this requires the simultaneous fulfillment of the two sets of equations 
([5]) and ([9|). We note that because of the dimensionality of A any solution of the previous 
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equations containing A and constructed with the only available (Lorentz-)covariant vector 
must necessarily be even under a change of sign of all coordinates x^ — >• —x^. Solutions odd 
in x^ exist but they require the involvement of parameters other than the coordinates and A 
(a wave vector, for instance, see section 4). 

The most general solution of this equation can be written as a superposition of both 
complex and real exponentials 

/d^ A/ / 71 \ 

^2 -j4 '^(^^ ~ (-^Mf ~^ ^f^i' ^™ ^ — ^ kx + B sinh kx)j — 7]^^, 

(33) 

with Ef^iy and D^y traceless, i.e. E^^ = D^^ = 0. In the previous expression E^j^^, Dfj,^, A 
and B are in principle all independent functions of k provided that the two following gauge 
conditions are met 

J (^^(^^ " ('^M^^^ sin kx + - J A sinh kx^ = (34) 

/ (04 ^^'^^ ~ ('^^^^ '^"'^ ''^ ~ '^^) " °' ^^^^ 

Clearly the integrands involved have to fall off sufficiently fast for large values of k for the 
integrals to exist. 

This solution has ten degrees of freedom to start with. Nine come from E^i, and D^i, after 
removal of the the trace. Another one comes from the coefficients A, B. Note that both A 
and B are needed to provide a full degree of freedom and likewise for E^^, and D^^. Using the 
gauge condition we can eliminate four of them, leaving six independent degrees of freedom. 
Unlike (I23p . the above solution does not admit any residual gauge transformation to further 
eliminate degrees of freedom. Any attempt to perform a residual gauge transformation would 
take the solution 'off-shell', i.e. the equations of motion would not be obeyed. 

On the other hand we have to ensure that h^i, « 1; However, in general this does 
not eliminate any degree of freedom, it is just a requirement of the linearized theory. This 
translates in requiring the first term in the expansion of the hyperbolic cosine to cancel the 
—r]^i, piece in (i33]l . or in other words 



/ 



"^^^ r6{k^ -2A)A{k) =4. (36) 



(27r)4 



Since (f33]) is the most general solution to the equations we must be able to recover the 
solutions in the previous section by performing an expansion in A. To do so we only have 
to choose the right form for E^y{k), D^y{k), A{k) and B{k). As mentioned previously, to 
reach a Lorentz-covariant formulation such as (j23p in the Lorenz gauge we can safely assume 
that D^y and B are zero as the resulting metric must satisfy h^y{x) = h^y{—x), as discussed. 
In addition A{k) can only be a constant on Lorentz covariance grounds. We will take it to 
be A{k) = ^ = Also E^^ needs to be a (traceless) Lorentz-covariant tensor, namely 
E^p{k) = ^ {kfikp — ^^A). The proportionality coefficient between E and A' comes from 
the gauge condition (j34p . Finally, as also indicated previously, the integrals require a finite 
support to be well defined and this should be implemented in a Lorentz-invariant way too; a 
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sharp cut-off will be used below, although this is not crucial at all. Expanding ([33]) . 

^^^^4 <^(fe^ - 2A) (^E^y{k) cos kx + ^A{k) cosh kx^ - rj^y 
= / - 2A) lyE„,{k) {l - *tll + . . (37) 



and using the definitions given above, 



r?^,^^ {k-x) 



(38) 



Now we introduce the cut-off, ^/2h.. Already condition (j36p dictates the value for A 

P 



/ _ 3277^ 



where C = -r L'^^ d\ k\ , . Then the solution reads 

■'U I ' , /r> A 1 



_ 2vr2 2y^ITI^ V 2A 2 7 2 4 AC 



2\/2A + A?2 

AC / ^ / A , 2 N A 27r2A 

=^ I " I 24 ^"^^""^ ^ ~ IQ^^"'' J ^ '^^'^'^ "C J ■ 

The value of E is fixed via the gauge condition ([31]) to E' = — -^f^, leaving the perturbation 
in the form ^ 

- (^M^f + '^'Hfj.ux'^) , (40) 

which is precisely (I23p . 



5 Wave-like solutions 

In this section we will finally investigate the effects of the cosmological constant in the prop- 
agation of GW in the appropriate coordinate system. 



5.1 Lowest order solutions 

We write /i^i/ = h^^ + /ij^. The term h^^ is the solution we just found, will be a 
perturbation on the metric induced by some source of GW. The same decomposition holds 
for the trace reversed metric h^y. Waves are usually considered in the transverse traceless 
gauge 

h^^ = h^'' = 0, = d^h^^' = 0. (41) 

This is compatible with the A-gauge condition as the r.h.s. of (jllh is unchanged when 
considering h^^ + provided that ([TT]) is fulfilled by h^^. This also makes clear that, at 
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this order, the gauge condition involves the perturbation associated to the background and 
not the metric perturbation associated to a gravitational wave. 

Since the proper equations of motion in the Lorenz gauge at this order, neglecting 
0{A.h^y), are just D/i^!/ = \^h^y + D^JX = 0) being the latter an independent perturba- 
tion, it is obvious that 



Uh 



w 



0, 



"M^^ - (42) 
and the gravitational wave solutions are in these coordinate systems functionally identical to 
those existing in flat space. 

Note that because the A/i^^ has been neglected, the remaining residual gauge invariance 
allows for a removal of four of the six degrees of freedom in /ij^ and the analogy with wave 
propagation in Minkowski space is complete. 

In the case of the lowest order equations the full solution of ()20p is 



h 



+ 2r]^yX^) + £')^ cos kx + D^^ sin kx 



(43) 



where = = 0, k^E^ = k^D'T = and k"^ = 0. 

We want to see now how plane waves such as the ones in (|43p look like in the new 
coordinate system. Transformation (jl7p acts both on the polarization tensors and on the 
arguments of the sine and cosine. For the polarization tensors we can always cut the expansion 
in A and keep terms only up to a certain order. However, the transformation on the arguments 
yields terms of the type Z^wA. which in general can be relevant. The sine and cosine can not 
be expanded, we have to transform the argument exactly; we shall later evaluate the error 
caused by retaining only the lowest order terms in the arguments. 

For the polarization tensors, since we transform them independently of the arguments, 
it is easy to see qualitatively what the corrections to the polarization tensors will be. On 
dimensional grounds alone, all corrections will be of order 0{^fA.Z) or at most ©(AZ^), being 
these quantities in the region of validity of the approximation very small. 

Nonetheless, the transformed wave-like solution to order \/A is 

/o 




(44) 



+ 0(A) 



The term w{T — Z) dominates the argument of the trigonometric functions and it can be 
checked numerically that the error made by omitting terms of order A or higher is ^ 10^^ 
for the purposes of next section. 
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5.2 Next-order solutions 



As we have argued before, it is not justified to neglect the term of order A/iJ^ in this case, as 
unhke for the case of the background, the magnitude of the two quantities is unrelated. We 
can add a wave-like piece to the solution (p3]) 



h 



(27r)4 



2A) ( E^y cos kx + D^u sin kx + {A cosh kx + B sinh kx) j — 7]^^ 



+ Ejf^ cos kx + sin kx. 



This will always be a solution of ([5]) and (j9]) as long as = D 



0, k.E'^^^ 



kuD u 



(45) 




and k'^ = 2A. However, now we are not allowed to perform any gauge transformation, at 
least at the next-order level. We can still use the gauge condition and the traceless condition 
to eliminate five degrees of freedom from the wave. We are left with a massive wave with 
five degrees of freedom. The polarization vectors of which, for a wave propagating in the z 
direction {ki = ^2 = 0), can be written as 



E, 



w 



w 

Vw'^-2A 
w ' 



E, 









w 


13 


Ell 


23 


El2 


00 


El3 



E- 



13 



Vw^-2A 



E' 



-Eu 



El2 
- Eqq-^ 

E23 



23 



2A 



Vw'^-2A 
El3 



Eoo\ 



-2 A 



E23 

„2 



-2A 



Ei 



(46) 



00 / 



And a similar expression for D^. At the exact level this is as far as one can go but in order 
to understand the meaning of these massive waves we turn again to an expansion in powers 
of A. We will proceed in two steps. First we expand the solution in powers of A and collect 
terms order by order. Then, using the same reasoning in the equations of motion, we can 
use an approximate residual invariance to rewrite the polarization tensors as the usual GW 
in Minkowski space-time plus an order A contribution with the extra degrees of freedom. 
The polarization vectors (j46p can then be written as 



E 



w 



/ Eqo 


El3 


-£'23 


-£'00 \ 


El3 


Ell 


-£'12 


-£'13 


E23 


E12 


—-£'11 


-£'23 


\Eoo 


El3 


-£"23 


-£^00 / 




El^J H 







+ 



13 





-Ae 

^E 
23 



y^2 El3 


ui-^ -£'23 


^■^oo\ 














P 2A 











2A rp , 



+ 0(A2) 



(47) 

The same decomposition applies to -D^. This expansion makes explicit the contributions of 
A at a given order. We want to expand 



(48) 



where the superscript refers to the order in A. The functions sine and cosine can also be 
expanded around a massless wave with coordinate-dependent amplitudes [15] 



h 



w 



E^ cos kx + sin kx 



E 



w 



^D]A coswit -z)+ + ^E^ ] sinwit 



(49) 
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or what is tantamount 



''fJ.U 



+ 0(A2) 



(50) 



We see that the massive wave we started with can be written at hnear order in the cosmological 
constant in terms of a massless wave where ah dependency in A appears only through the 
polarization tensors 



h 



w 



cos w^t 



(51) 



where and can be read from 



50]). The above is a valid solution of D/i^ + 2A/i^ = 
only to order A (included), which means we can expand the equations of motion to the same 
order without loss of validity 



□/iJS + n/iW + 2A/iW + 0(A2) = 



(52) 



Now we can split the problem and solve order by order 



□/.W = 
□/i«+2A/ig)=0 



(53) 



Due to the fact that ()52p is not exact, the solution to it can admit a residual gauge transfor- 
mation that will take the solution 'off-shell' some order beyond the order we consider. For 
the transformed solution 

Uh'IS) = 



The first equation in (j54p is analogous to (j42p . i.e. residual transformations on h^^} are not 
restricted. To order zero we obtain GW analogous to the ones in fiat space (in the present 
set of coordinates, that is). But in this case the transformation propagates to the following 
order through the second equation in (154p making necessary to find the transformed h!l]^ . 
It is not difficult to see that the following polarization tensor fulfills the necessary require- 



ments of tracelessness as well as the gauge condition {k^E 



0) 



/ 7^E( 



E 



w 



00 



V 



tE 



13 



-E'23 

Eqq 



Eu - 

E12 — -7^Di2 



11 



w^E23 

Ei2 - 



A 



-E 



11 



12 

+ —^11 

I -7/1 -L-L 



A 
— ^ 



E 



13 



A 
— ^ 



E. 



23 



j.Eqq \ 

^^13 
7E23 
Eqo 



(55) 



Dfj^u is similarly obtained from (|50p . Notice the presence of the usual components (of 0(1)) 
in the polarization tensor in the x, y entries of the metric. 

To this order in A we obtain massless waves with coordinate-dependent modified ampli- 
tudes which depend on A. We can see that the extra degrees of freedom due to the form of the 
linearized equations of motion for non-zero A will only couple to matter fields proportionally 
to A thanks to the coupling li^^T^^ and thus will be irrelevant in practice. 
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5.3 Transformed next-order solutions 



Now we are ready to apply the series of coordinate transformations l\27\ [29l [311 [T7|) to the 
wave-hke solution (I5ip that we found in the previous subsection in order to obtain a physical 
expression in FRW coordinates. Recall the waves in the general Lorenz gauge read 



h 



w 



E^^{A,z) cos w{t - z) + D^^{A,z) sin w{t - z), 



w I 



(56) 



where can be read off from (|55p . From (|56p it is clear the only modification with respect 
to the plane waves of the lower order is in the polarization tensors, being already of order A. 
This suggests that all the new modifications to order A of the next-order waves are due to 
the change of coordinates. Explicitly the transformed waves to order A read 






o\ 


5A 72 \ 











+ 0(A3/2) 



^wA (T^ + T'^Z - 5TZ^ + 2Z^) + ©(A^/^)^ 



At-1 I 2At^2 I 5A 72 
3i + — i + jgZ^ 



I\rp I Zl\rp2 I 5A r72 
-ji + -g-i + JgA 






0/ 



+ 0(A3/2) 



^wA (T^ + T^Z - 5TZ^ + 2Z'^) + ©(A^/^)^ . 

(57) 



6 Detectability 

Let us now do some order-of magnitude estimates to evaluate the effect of the corrections 
induced by A 7^ on the propagation of gravitational waves. 

For the polarization tensors we have not attempted to derive the A-order corrections in full 
detail, although this is possible, because already the most relevant correction, i.e. \^ZE^J , 
has to be some orders of magnitude smaller than Ej^J for the approximation to be valid. For 
example for a coordinate value of the order of a typical distance to a supernova, 1023 m, the 
quantity VAZ ~ 10-3 (A ~ 10"^2 ^ g-2^)^ rj,^jg already means a small correction 
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to an amplitude that has so far escaped detection and which presumably will not be measured 
with sufficient precision to discern the effect of the A-order effects in the foreseable future. 
However, conceptually it is an interesting result. 

It is more interesting to work out the corrections to the dispersion relation for (I44p . As 
previously, let us consider waves that propagate in the Z direction and are monochromatic. 
The maxima of the wave will be reached when 



TZ 



nvr, 



or 



Z^i,^{n,T) = T 



nvr 
w 



rp2 



A n^TT^ 
3 



From (I59p we can also calculate the phase velocity of the wave which is defined as 



VpiT) 



dZr, 



dT 



1-T\I- + 0{A). 



(58) 
(59) 

(60) 



We see that in comoving coordinates the phase velocity is smaller than 1. This does not mean 
that the waves slow down. We can calculate the velocity in 'ruler' distance. For a fixed time 
we have 

I A 

\- _ (61) 
I A . .„ ^ 

It is also interesting to rewrite the trigonometric functions of the wave defining Wef[{Z) 
w ( 1-Z./4 




cos 




COS W(.q{T — Z) 



(62) 



Note that the transformed wave corresponds to a usual wave with an effective frequency 
dependent on the coordinate Z. The wave becomes red-shifted as it propagates away from 
the source. 

To see explicitly the effect of A in the propagation of a wave described in comoving 
coordinates we plot (Figure 1) one of the components of the wave for a given instant 
(r = for simplicity). A wave with a physical frequency ranging lO^Hz < w < 10^^'^Hz 
cannot be practically plotted in the relevant Z-range. To see the effect in a few cycles we 
take w = A - 10~^^Hz, which does not affect the overall magnitude of the correction. We plot 



the wave for A = 10" 



-52 



m 



and for A = 10" 



^^m ^ to assess the influence of A on the wave 



-Zw 1 



7. A 

^a/ 3 



propagation. Then we plot ~ (l + |AZ^) cos 

From these results we can already draw some conclusions. The genuine corrections due to 
the mass-like term in ^ remain unchanged in the transformed waves if we cut the expansion 
to order 0{A). Moreover they are of order which is in practice irrelevant unless the 
value of A is much greater than the current value. However, transformation (|17p induces 
modifications to the wave, both in the amplitude and the phase, of order \/A and A. This 
modifications result in a simulatneous increase of the wave-length and of the amplitude with 
the coordinate Z. As shown in Figure 1, the most interesting region for detection would be 
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Figure 1: Dependency of the amplitude and wave-length on the coordinate distance Z (ex- 
pressed in meters) for a constant value of T and for different values of A: The dashed line 
corresponds to A = 0, the dotted line to A = 10~^^m~^ and the solid line to A = 10~^^m~^. 



that of events (supernovae and pulsars for example) happening at a distance Z ~ 10'^'^ — lO'^^m 
away, for which the correction \J\z ~ 10~^ — 10~^ is not negligible and is well within the 
validity range of the approximation. In fact to have this type of correction into account seems 
probably essential to properly account for the measurements of this type of phenomena in 
pulsar arrays. 



7 Summary 

The purpose of this work was to investigate the effect of the cosmological constant in the 
propagation of gravitational waves in a linearized theory of Gravity. The presence of A leads 
unavoidably to the curvature of the background space-time in which the waves propagate. 
Within the linearized approximation (which is the only framework where one can properly 
speak of 'waves') this leads to a decomposition g^^y ~ + + including a modification 
of the background (corresponding to the curvature) and a wave-like perturbation. 

To see the way the propagation of the waves is affected, one has first to understand the 
implications that the different coordinate choices (gauge choices) have in the resolution of the 
equations of motion as well as the importance of the terms of different order retained in the 
linearization. One is free to choose any particular gauge to solve the equations, however since 
the linearized Einstein equations are not invariant under general coordinate transformations 
their form will depend on the gauge choice. We argue that the above procedure of linearization 
is consistent in some coordinate systems but not in others. In particular, it is inconsistent to 
linearize the equations in the familiar Priedmann-Robertson- Walker cosmological coordinates 
(the metric only depends on time). 

Einstein equations can however be consistently linearized in Schwarzschild-de Sitter co- 
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ordinates; then h^^ corresponds to a linearized version of the SdS metric, expanded to first 
order in A. This metric can be easily modified to fulfill the Lorenz gauge condition. In this 
particular gauge, i.e. in this particular choice of coordinates, the analysis of gravitational 
waves follows a pattern very similar to the one in Minkowski space-time. In the case where the 
A/i^jy term is dropped the residual gauge freedom of the Lorenz gauge allows for the removal 
of four additional degrees of freedom in the general solution, leaving the wave-like component 
with the usual two physical degrees of freedom of waves propagating in flat space-time. 

On the contrary, if the term Ah^u is retained in the equations of motion the situation 
changes. Even in the Lorenz gauge the invariance under residual gauge transformations 
is lost. Again it is not hard to find the most general solution to the linearized equations 
composed of a background and a wave-like components. We prove the background solution 
to be consistent with the result previously found if A is small. Since there is no residual 
invariance, the wave-like solution has to be interpreted as a massive wave with five degrees of 
freedom (the gauge condition and the trace condition amount to five constraints). However, 
we can make use of the approximate residual invariance at the leading order in A to rewrite 
the solution as massless gravitational waves with position-dependent modified amplitudes 
that change very slowly given the current values of A. There are only two 0(1) polarizations; 
the remaining degrees of freedom (up to the five independent ones required for a massive spin 
two wave) are of 0{A) and couple extremely weakly to matter sources. 

Finally, one has to transform these solutions to the physically significant FRW coordinates 
in order to extract observable consequences. At this point modifications of 0{^/A) appear. 
Numerically these can be quite relevant for certain gravitational waves travelling from far 
away sources and the effect of A can absolutely have a detectable impact on pulsar timing 
arrays. Waves are modified both in the phase and the amplitude; in cosmological coordinates 
they are red-shifted in a prescribed way and the amphtude of plane waves grows as they move 
away from the source. 
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